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EFFECTS OF CONCENTRATION AND VIBRATIONAL RELAXATION 


ON INDUCTION PERIOD OF HYDROGEN- OXYGEN REACTION 
by Frank E. Belles and Milton R. Lauver 
Lewis Research Center 


SUMMARY 

The possibility that slow vibrational relaxation of oxygen (Og) or hydro- 
gen (Hg) might affect the induction period of the high-temperature Hg - Og re- 
action is considered in terms of the data from three different kinds of shock- 
tube experiments. The data cover a wide range of initial composition and tem- 
perature and include observations in Hg - air,, Hg - Og - argon (Ar), and Hg - Og 
mixtures. It is found that induction times in excellent agreement with experi- 
mental data can be calculated by solving the differential rate equations based 
on the following simple kinetic scheme: an initiation reaction between molec- 

ular species, plus the three usual chain-branching reactions. It is not neces- 
sary to invoke vibrational effects to obtain this agreement, but it is neces- 
sary to solve the problem without incorporating steady-state assumptions or 
limitation on relative concentrations of Hg and Og . Thus, it is concluded that 
existing data exhibit little or no effect of slow vibrational relaxation. In 
support of this conclusion, calculations are presented to show that the greatest 
possible effect of slow Og relaxation would, in fact, be a small one. The cal- 
culations were carried out by making the extreme assumption that all the activa- 
tion energy of the slow chain-branching step must be supplied by vibrational 
energy of Og. Since the level of this energy varies during the induction period 
in a manner governed by the vibrational relaxation time, the reaction rate con- 
stant becomes time dependent. Unfortunately, the equivalent calculations for Hg 
cannot be made because the relaxation time of Hg in the presence of Og is not 
known. 


INTRODUCTION 

If a gaseous mixture containing hydrogen and oxygen is suddenly heated to 
a temperature greater than about 1000° K at moderate pressure, for example, by 
means of a shock wave, its ignition is preceded by a short induction period. 
During this interval chain branching occurs, and since the induction period 
comes before the heat -releasing part of the reaction, the chain-branching pro- 
cess takes place under conditions of essentially constant temperature and pres- 
sure. Moreover, the time is so short, a few to a few hundred microseconds, 
that wall effects cannot make themselves felt. 



Thus, measurements of the induction period permit the chain -branching pro- 
cess to be studied in a relatively uncomplicated environment. This fact was 
perhaps first recognized by Schott and Kinsey (ref. 1), who made an extensive 
set of measurements over a wide range of temperat lores and values of the hydro- 
gen to oxygen concentration ratio [H 2 ]/[0 2 ] from 0.25 to 5. 

Only four reactions are needed to describe the buildup of free radicals 
during the induction period. First, there must be an initiation reaction in- 
volving molecular species. In the lower range of temperatures (1000° - 2000° K) 
this reaction may be the following (ref. 2): 

h 2 + Og 2 OH (i) 

After a short initiation period, which calculations show will occupy only a 
small part of the total induction period (ref. 2), the following reactions 
branch the chain: 


OH + H 2 - H 2 0 + H (I) 

H + 0 2 - OH + 0 (II) 

0 + H 2 -> OH + H (III) 

Rate constants k* , kn , k 2 , and k* are available for all of these reactions 
(refs. 2 and 3) 1 1 

Schott and Kinsey correlated their data by plotting the parameter t^tOg] 
against reciprocal temperature. They took the induction period tj,, -which of 
course must always be defined somewhat arbitrarily, as the time at which the 
hydroxyl concentration [OH] reached the limit of detectability in their experi- 
ments, about 10~^ mole per liter. They noted that this correlation is the one 
to be expected if chain branching is controlled by the rate of the slow step 
(reaction (il)). 

However, they observed a trend in their data such that tj_[02>] increased 
with increasing mole fraction of oxygen at constant temperature. This is not 
anticipated on the basis of the simplified kinetic model in which the control- 
ling rate is that of reaction (il). Schott and Kinsey tentatively ascribed 
this trend to slow vibrational relaxation of oxygen, which might tend to reduce 
the rate of reaction (il) below the value it would have under conditions of 
full thermal equilibrium. This amounts to the suggestion that at least part of 
the activation energy normally resides in vibrations of the oxygen molecule. 

In order to investigate this important possibility more fully, it is nec- 
essary to consider data for the widest possible range of [ H 2 ] / [ 0g ] values and, 
hence, of relaxation times; the data of White and Moore (ref. 4) are extremely 
valuable on this account. If the data are to be interpreted in the light of 
reactions (i) to (ill), then it is also important not to suppress any implica- 
tions of this kinetic scheme; that is, the differential rate equations must be 
solved without incorporating any assumptions about steady states or limitations 
on [ H 2 ] / [ O 2 ] ratios. If, after this, any discrepancies remain, they can be 
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terms of slow vibrational relaxation. 


EFFECTS OF CONCENTRATION ON INDUCTION PERIOD 
Calculat ions 

In order to calculate the concentration of H, 0, or OH as a function of 
time during the induction period, it is necessary to integrate the set of dif- 
ferential equations based on reactions (i) to (ill). An obvious way to solve 
the problem without making any simplifying assumptions is by numerical inte- 
gration carried out with the aid of a high-speed computing machine. Such nu- 
merical integrations have been reported (ref. 5) for 5-percent-Hg - 95-percent- 
air mixture, and they were used to explain the chemiluminescence of hydroxyl 
radical during the induction period. These calculations showed that after a 
brief initiation period the free-radical concentrations adopt the following 
simple time dependence: 


EL = EL = ML = exp (t /t ) 

r tt -i r /"\ n ta tt t ' 9 


[H] 


[0] Q [OH] 


(1) 


where t is the time and t is the exponential time constant. The quantities 
with subscript 0 are pseudoinitial concentrations obtained by extrapolating 
the linear part of a semilogarithmic plot to time zero. 


The simple form of equation (l) strongly suggests that there is a general 
analytic solution for the induct ion -period kinetics. Actually there is a 
general solution due to the fact that reactions (i) to (ill) are in effect 
first-order reactions because [0g 1 and [Hg ] remain constant during the induc- 
tion period. The method of solution for such situations is outlined in stan- 
dard text books (ref. 6), and, the solution was carried out by Kondratiev 
(ref. 7) for the kinetic scheme used in the present paper. Nicholls, Adamson, 
and Morrison (ref. 8) also appreciated the inherent simplicity of the 
induction-zone kinetics, but they dealt with a more extensive set of reactions 
and introduced approximations that lead to a solution good only if [Hg] and 
[0 2 ] are of order unity. Their solution, therefore, is not adequate to deal 
with a very large range of [Hg]/[0g]. 

Brokaw (ref. 9) has recently presented an approximate technique by which 
an analytic solution can be obtained that gives excellent agreement with the 
results of numerical integration of the induct ion- zone kinetics. Hfe obtains 
the following expressions, which are applicable to the high-temperature regime 
of short induction times: 


7 = [H * 3 | 


8 V° 2 ] 

k,[H P 


1/2 


- 1 


2 + 1 + 


v° 2 ] 


1/2 


8k ? [0 ? A 1 / 2 

1 + — = — —I + 1 




( 2 ) 


3 


and 


[0HI 0 


2 M°2 ] 


(3) 


Equation (3) expresses the reasonable notion that the initiation reaction will 
cease to be important when the fastest chain-branching reaction (i) uses up 
OH as rapidly as it is produced by the initiation reaction. 

From the previous results the induction time of reference 1 ([OH] = 10”^ 
mole/liter) multiplied by the oxygen concentration can immediately be written 
as 


[ ° 2 ] 

t, [0 ? ] = — — < 

1 2 rV 


L + k iM 

'll + 8 k 2 [ 0 2 ] \ l/2 + x 

\ k 2 ^° 2 V 

A k 3 [ W J 


1/2 


2k l k 3 


1 + 


3k 2 [0 2 ] 


, 1/2 


In 


10 


-6 




J 


(4) 


There are two main points to be made about equation (4). The first is 
that the parameter t [0g] should not succeed in correlating induction times 
obtained over a wide range of [H 2 ]/[0 2 ] values. The second is that the effect 
of the initiation reaction is small; since it is in the logarithmic term; thus, 
a large error could be made in [0 H] q without materially affecting the result 
calculated from equation (4). This is fortunate., because there is no indepen- 
dent evidence as to the rate constant or even the occurrence of reaction (i). 


We also want to compare calculated induction times with experimental 

values obtained by two other techniques. The most important set of data is 

due to White and Moore (ref. 4); who measured the period of constant density 

behind shock waves by interferometry. In a preliminary report (ref. 10), White 

noted that the times obtained in this way are similar to those reported in 

reference 1. It is assumed, therefore, that the two kinds of measurements are 

essentially equivalent; equation (4) is also used for comparisons with White 

and Moore's pesults. 

* 

Finally, a limited amount of new data will be presented, obtained from 
observations of light emitted at 3080 Angstroms by electronically excited 
hydroxyl radical OH*. The interpretation (ref. 5) of this light predicts an 
inflection point in the curve of light intensity against time behind the shock 
front because of the competition between the reaction producing the excited 
molecule 


H + 0 2 + H 2 - H 2 0 + OH* 


(W) 
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and collisional quenching by water produced in reaction (l), 


OH* + H 2 0 -> OH + H 2 0 (V) 

The inflection time is given by the following expression (ref. 5): 


t inflect ion T ln 

The quantities (i/t) and [OH]q may be obtained by means of equations (2) and 
(3). The rate constant k^ is taken as the binary collision number. 

The other rate constants used for calculations are listed in the following 
table where k = A exp(-E/RT). 


1 (i Zilf 


l 


(5) 


Reaction 

Preexponential factor. 

Activation energy. 

Source of 


A, 

E, 

data 


liter/mole-sec 

cal/mole 


i 

1.QX10 11 

70,000 

Ref. 2 

I 

6.3X10 10 

5,900 

Ref. 3 

II 

4.QX10 11 

17,000 

Adjusted 0, 

III 

1.2X10 10 

8,950 

Adjusted 3, 


a Adjusted to give good fit to light -emission data obtained for 
5-percent -H 2 - 95-percent -air mixture over the temperature 
range from about 1200° to 1600° K. These k's represent 
increases in the values given by Baldwin (ref. 3) for k£ 
and by Kaufman and Del Greco (ref. 3) for k^ . The maximum 
change is a factor of 1.8 times Baldwin's k 2 

(2.0X10 11 exp( -16600/RT) ) and 3.4 times Kaufman and Del Greco's 
k 3 (2.5X109 exp( -7700/RT) ) . 


They are either literature values or are adjusted with literature values as a 
starting point. In effect, all comparisons between calculated and observed 
induction times have been normalized to one particular set of experimental 
data by means of the tabulated set of rate constants. However, it should not 
be inferred that this set is unique or that it represents improved absolute 
values. 


Comparison with Experiment 

Light -emission experiments . - Inflection times were measured from oscillo- 
scope records of light intensity behind shock waves traveling through 5- and 
20-percent H 2 - air mixtures at 10 torr Initial pressure. Zero time was fixed 
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Inflection time, Msec 


by the arrival of the shock at a thin-film resistance gage located at the same 
axial position as the slit through which emitted light was transmitted to a 
monochromator and photomultiplier detector. Details of the experiment are 
given in reference 5, which also includes the data for the 5-percent mixture. 
Results are plotted in figure 1, which also shows inflection times calculated 

from equation (5) as solid 



lines. The temperatures 
used to plot the data and 
the solid lines were cal- 
culated from shock veloc- 
ity by using the graphical 
method of Markstein (ref. 
11), and they represent 
full thermal equilibrium 
behind the shock. 

Since the rate con- 
stants were chosen for a 
good fit to the 5-percent 
data at intermediate tem- 
peratures, the upper solid 
line naturally agrees well 
with the experimental re- 
sults. Inflection times 
for the 20-percent mixture 
are predicted to be about 
half as great at a given 
temperature. Despite a 
higher degree of scatter, 
due largely to the fact 
that shock velocities were 
not as uniform as in the 
leaner mixture, the 20- 
percent data bear out this 
prediction quite well. 

It was pointed out in 
the earlier work (ref. 5) 
that equation (5) is de- 
rived from an approxima- 
tion that becomes worse as 


Temperature, T, °K 


the temperature increases. 


Figure 1. - Comparison of calculated and observed times of inflection in records of 
OH* 2 I + - 2 n emission intensity as function of time. 


In order to assess this 
effect, the differential 
equations based on re- 
actions (i) to (V) were 


integrated numerically to give [OH*] as a function of time, and inflection 


times were obtained from the plotted results. These more accurate times were 
only 10 percent larger at the highest temperature considered (1900° K), so 
equation (5) is indeed a very good approximation. 
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The simplifying kinetic assumptions adopted by Schott and Kinsey (ref. l) 
and by Nicholls, Adamson, and Morrison (ref. 8) lead to the conclusion that 
(l/O should be given by the following expression: 

1/t — 2k^ [Og ] (6) 

This result, when inserted in equation (l), provides the basis on which Schott 
and Kinsey expected their data to correlate in terms of the parameter t^tOg]. 
It can be seen from equation (5) that the same correlation should hold for the 
light-emission data of figure 1 if equation (6) were correct; but [Og] is only 
about 1.2 times as large in the 5-percent as in the 20-percent mixture, while 
the data differ by a larger factor; so the correlation will not bring the re- 
sults together. 

Light -absorption experiments . - Although the work of Schott and Kinsey 
covered a much wider range of [ r Hg]/[Og] values than the light -emission studies 
just discussed, the failure of the tjJO;?] correlation was obscured by scatter 
in their results. However, they did observe that when the data for each mix- 
ture was considered separately, t_^[Og] varied with mole fraction of oxygen. 

The values they obtained for a constant induct ion- zone temperature of 1800° K 
are quoted in the last column of the following table: 


1 

2 

3 

4 

5 
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Mix- 

ture 

°2 > 

[h 2 ]/[o 2 ] 

a [0 2 ], 

( mole/liter )X10^ 

V 

secXlO^ 

t i 

[ 0g ] , sec ( mole/liter )xl0^ 

percent 



Calculated 

Observed 






Constant kg 

Time- 
dependent , 

k 2 







E<1. (4) 

Numerical 


A 

0.43 

1.6 

2.5 

16.2 

5.7 

5.7 

6.1 

2.7 

B 

.49 

2.0 

2.6 

15.1 

5.1 

5.4 

5.7 

3.7 

D 

.45 

5.5 

4.0 

4.6 

3.3 

3.6 

3.8 

3.1 

F 

.49 

2.0 

12.0 

3.2 

4.5 

4.6 

4.9 

4.5 

C 

1.99 

1.9 

12.0 

4.0 

4.6 

4.7 

5.2 

4.7 

E 

2.00 

.5 

21.0 

7.3 

8.9 

9.2 

10.7 

5.4 

H 1 

4.00 

.25 

37.0 

8.1 

13.5 

13.8 

16.5 

7.8 

H 

19.7 

.25 

49.0 

7.7 

13.2 

13.7 

17.4 

11.7 


Unpublished data received from G. L. Schott, Dec. 1963. 

Each mixture is designated by a letter as in the original work (ref. l). 

The first two columns of calculated results give values obtained by means 
of equation (4) and by numerical integration. There is no more than 10-percent 
difference in any of the pairs of numbers, so it can be concluded that equa- 
tion (4) is valid over the full range of [HoJ/tOg] values. Comparison of 
either column of calculations (constant kg) with the experimental results 
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Temperature, T, °K 


<b) Richer mixtures. 

Figure 2. - Comparison of calculated induction times with experimental 
data of White and Moore. 


is very gratifying. The worst 
discrepancy is a factor of about 
2 for mixture A. Thus, simple 
calculations based on reactions 
(i) to (ill) can explain the re- 
sults of the light-absorption 
experiments quite well without 
invoking effects due to slow vi- 
brational relaxation. 

Interferometric experi- 
ments . - The induction-period 
measurements of reference 4 
cover an extremely wide range of 
[Hgj/tOg] values (0.0075 to 24), 
so they provide the best test of 
the simple kinetic scheme adopt- 
ed in the present paper. In 
figure 2 the data are reproduced, 
together with lines calculated 
for eight of the nine [h^l/tOg] 
values studied. The calcula- 
tions were made by means of 
equation (4) by assuming a con- 
stant induct ion -zone pressure 
of 1 atmosphere in all cases. 
Inasmuch as the pressure only 
enters via [Og], in the loga- 
rithmic term of equation (4), 
the result is quite insensitive 
to the assumed pressure, and 
moreover (as shown by an un- 
published communication from 
White) 1 atmosphere is a reason- 
able average for these experi- 
ments. 

Inspection of figure 2 
shows that experiment and cal- 
culation agree very well indeed. 
At a given temperature, the 
parameter t^tOg] decreases as 
the ratio [Hg]/[0g] increases, 
spanning nearly two orders of 
magnitude from the leanest to 
the richest mixture. The major 


discrepancy seems to be in the 

richest mixtures, where calculations for [H 2 ]/[0 2 ] values of 7 and 24 give two 
lines quite close together; whereas, the data for these two mixtures appear to 


be more widely separated. It must be borne in mind, however, that the calcu- 
lated tj_ corresponds to the time at which [OH] = 10“ 6 mole per liter, and 


8 



there is no real assurance that the data of White and Moore should coincide, 
since they are based on density changes rather than on the concentration of OH. 
In any event, the discrepancies between experiment and calculation are no more 
than a factor of two at worst. 


The slopes of the calculated lines undergo a noticeable increase as 
[h 2 ]/[o 2 ] increases, a rather marked change taking place between [HrJ/tOg] 
values of 0.10 and 0.33. This behavior is due to the increased importance of 
reaction (ill) with its smaller activation energy in the leaner mixtures, while 
in those mixtures less abundantly supplied with 0g, reaction (il) assumes more 
control. 


White and Moore (ref. 4) observed that their data can be correlated if the 
parameter t j_( [0-> 3 [Hg ] is used instead of t^tOg]- It is not readily ap- 

parent that equation (4) can predict this correlation. It does follow, how- 
ever, if one considers a less exact solution of the induct ion -zone kinetics. 
Brokaw (ref. 9) found that a rather good first approximation to (i/t) can be 
obtained by imposing a steady state on OH (but not on 0 and H) : 


1 

T 


kjlHj] 

2 



8 k 2 [ 0 zA l/Z 

/ 


l 


( 7 ) 


The approximate time at which [OH] = 10“^ mole per liter is therefore 



/ + 8k 2 t0 2 ]\ l/2 ^ 

-1 

“ — 

c 2 

In 

10" 6 

k 3 [H 2 ] 

\ ) 

_2k i [0 2 ]/k 1 _ 


( 8 ) 


Equation (8) shows that the parameter 
correlate the data provided that 


t i ([0 2 ][H 2 l) 1 / 2 


will approximately 


(a) In 


10 


-6 


2k i [0 2 ]/k 1 

1/2 


— constant 


, v / 8k [0 ] 

Ob) 1 + £_ 

V k 3 ^ H 2 ^ 


» 1 


At any given temperature in the experimental range, calculation shows that 
criterion (a) is well satisfied for all [Hg]/[0g] values studied by White and 
Moore and that criterion (b) holds very well for the leaner mixtures up to 

[Hgl/tOg] = 0.33. Consequently, a plot of t^( [ 0g] [Hg] )^ 2 as a function of 
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l/T should, at the very least; “bring together the data for all the lean mix- 
tures; hence; the correlation does have a basis in the induct ion -zone kinetics. 


DISCUSSION 

The foregoing comparisons between calculated and observed induction times 
show that the experimental results are explained by straightforward application 
of a simple kinetic scheme. There is no need to invoke effects due to slow 
vibrational relaxation of 0 % or H->. 

Perhaps this is not a surprising result; because there is no a priori 
reason to believe that reactions (i) to (ill) normally obtain part of their 
activation energy from internal degrees of freedom. On the other hand; many 
reactions do. The measured rates of thermal dissociation of simple molecules; 
for example; usually have to be interpreted by assuming that internal energy is 
a contributing factor (ref. 12). Of course; such dissociation reactions are 
strikingly different from those occurring in the induction period of the 
Hg - Og reaction; but it is nevertheless worthwhile to see if calculations can 
confirm that vibrational effects should be absent; or at least small. 

First of all; it must be acknowledged that some or all of the initiation 
process; whether by reaction (i) or some other reaction; will proceed in vibra- 
tionally cold gas. But equations (l) and (5) show that induction time is ex- 
tremely insensitive to the rate of initiation; therefore; it is unlikely that 
experimental investigation could detect vibrational effects -here. 

Of the much more important reactions; namely; those that branch the chain; 
only reaction (il) involves the oxygen molecule. Therefore; let us next con- 
sider what would happen to the induction time if the vibrational energy of 0 
were needed for activation. The simplest way to compute an upper limit for the 
effect is to assume that all the activation energy must reside in vibrations. 

In other words; RT is replaced by RT y ; where T is the overall temperature 
of the induction zone and T v is the vibrational temperature of Og at a given 
instant; yielding 


k 2 = A 2 exp(- E 2 /RT v ) (9) 

The timewise history of starting from 300° K, may he written as follows 

(ref. 13): 


T y = T - (T - 300 ) exp ( - t/t y ) (10) 

where t v is the vibrational relaxation time. 

White and Millikan have presented data (refs. 14 to 16) and a general cor- 
relation (ref. 17) that permit the relaxation time of O 2 in any mixture con- 
taining O 2 , H2j Argon (Ar), and nitrogen (N 2 ) to be estimated. Values calcu- 
lated at 1800° K for the H 2 - O 2 - Ar mixtures used in reference 1 are listed 
in column 5 of the table on page 7. These times differ considerably from those 
estimated by Schott and Kinsey because of the subsequent discovery (ref. 16) 
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that Eg is extremely effective in collisionally exciting the Og vibrations. 

The range of ratios of observed tj_ to t v is 2.6 (mixture H f ) to 16.9 (mix- 
ture D); hence, the time required for relaxation is always less than the mea- 
sured induction times. 

The differential rate equations containing the time -dependent kg 
(eqs. (9) and (10)) were integrated by means of a high-speed machine program. 
The time at which [OH] = 10 ”6 m ole per liter was read off from the plotted re- 


sults. Figure 3 contrasts the behavior 



Figure! - Effect of time -dependent k 2 on hydroxyl concen- 
tration during induction period in 0.01 hydrogen - 0.04 
oxygen - 0.95 argon mixture at 1800° K. 


behind the shock will be about 2070° K. 
Schott and Kinsey were (except for H) 
ess of vibration relaxation had little 
temperature . 


when kg is constant and when it is a 
time-varying function of the Og vi- 
brational temperature in mixture E ! , 
the mixture for which t v comes 
closest to t^. Remarkably, the ex- 
treme assumption that all the activa- 
tion energy must come from the vi- 
brational energy of Og only increases 
tj_ by 20 percent. Similarly small 
effects for the other mixtures are 
shown in the table on page 7. 

The slightly more difficult 
integration based on reactions (i) to 
(V) was also carried out subject to a. 
time-varying kg so that inflection 
times could be obtained for compari- 
son with the lighi; -emission experi- 
ments in 5- and 20-percent Eg - air 
mixtures. The results are plotted 
as dashed lines in figure 1, page 6. 

At this fine it is necessary to 
decide whether the data points in 
figure 1 can be compared with the 
dashed lines . It will be recalled 
that the points represent observed 
inflection times plotted against the 
reciprocal of the full- thermal- 
equilibrium temperature, as calcu- 
lated from shock speed by Markstein' s 
method (ref. ll). But there is actu- 
ally a marked temperature variation 
in the induction zone, because these 
mixtures consist of diatomic gases 
that absorb energy in vibrations. 

For example, if the thermal- 
equilibrium temperature of the induc- 
tion zone in 20-percent Hg - air mix- 
ture is 1900° K, the temperature just 
In contrast, the mixtures used by 
heavily diluted with argon, so the proc- 
effect on the overall induct ion -zone 
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However, calculation of the vibrational relaxation times for 0 2 in the 
H 2 - air mixtures shows that the process is very fast, because of the appreci- 
able mole fraction of H 2 present. At 1500° K, for instance, the observed in- 
flection times are about 6 times longer than the calculated relaxation times in 
the 5-percent mixture, and in the 20-percent mixture, about 10 times longer. 
Corresponding calculations for the N 2 present in the mixtures suggest that 
relaxation will require times comparable to the observed inflection times; in 
actual fact, however, the process will be greatly accelerated (ref. 4) by ex- 
change of vibrational energy from the easily excited 0 2 . Therefore, the tem- 
perature of the induction zone during most of its history will be close to the 
full- thermal- equilibrium value, and it is indeed appropriate to use that tem- 
perature in plotting the data. 

In figure 1, page 6, it is seen that the points for higher temperatures 
perhaps come closer to the dashed lines calculated for the case of slow C>2 
relaxation. However, the predicted effects are again small, just as in the 
mixtures studied by Schott and Kinsey, and the data are certainly not good 
enough to warrant any positive conclusion. 

It is more difficult to consider what might happen if slow vibrational 
relaxation of H 2 played a role in the chain-branching kinetics. In the first 
place, H 2 participates In two of the reactions (i and III), and it is best to 
consider only one at a time for the sake of clarity; but more Important is the 
fact that the vibrational relaxation time of H 2 in H 2 - 0 2 mixtures is unknown. 

If relaxation is assumed to occur by a simple collision process, the time 
can be calculated (ref. 17), and it turns out to be extraordinarily long due 
to the very high characteristic temperature of the molecule. In fact, the 
collisional relaxation times turn out to be much longer than observed induc- 
tion periods. However, White and Moore (ref. 4) have found that the process 
is actually a great deal faster than it is expected to be on the basis of the 
correlated data for other systems (ref. 17); in fact, it is so fast that it 
cannot be resolved interferometrically in their experiments. Therefore, values 
of ty for H 2 in the presence of O 2 are unavailable at present, and it is 
not possible to carry out calculations in which equation (10) is used to make 
kq or k 3 time dependent. But reasoning by simple analogy from the results 
obtained with a time -dependent k2, it certainly seems likely that the effect 
of H 2 relaxation on induction time would be small. 


CONCLUSIONS 

From this study on the effects of concentration and vibrational relaxation 
on the induction period of hydrogen (H 2 ) - oxygen (O 2 ) reactions, the following 
conclusions are drawn: 

1. Induction times for the hydrogen- oxygen reaction, measured behind 
shock waves by various techniques and over a very wide range of [H 2 ]/[0 2 ], 
can be explained by a simple kinetic scheme consisting of an initiation re- 
action and the usual three chain-branching reactions. 
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gen (N 2 ) does not play a significant part in the induction-period kinetics. 
While perhaps not unexpected, this fact is pertinent to practical considera- 
tions of air-hreathing hydrogen-fueled engines., , , - / 

/k/1^ r e\A A te^ . 

4. Calculations based on extreme assumptions about the possible role of 
vibrational energy in the chain-branching reactions show that slow relaxation 
of Og could have only a small effect on induction time, probably not detectable 
experimentally. Equally confident statements cannot be made about the possible 
results of slow relaxation of Hg because its relaxation time in H 2 - 0 2 mix- 
tures is unknown; however, it is very likely that this effect (if it exists) is 
also a small one. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, May 21, 1964 
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